arXiv: 1505.05913vl [math.DS] 21 May 2015 


Dynamics of a generalized Beverton-Holt competition model subject to 

Allee effects 

Yun Kang 1 


Abstract 

In this article, we propose and study a generalized Beverton-Holt competition model subject to Allee 
effects to obtain insights on how the interplay of Allee effects and contest competition affects the persis¬ 
tence and the extinction of two competing species. By using the theory of monotone dynamics and the 
properties of critical curves for non-invertible maps, our analysis shows that our model has relatively 
simple dynamics, i.e., almost every trajectory converges to a locally asymptotically stable equilibrium 
if the intensity of intra-specific competition intensity exceeds that of inter-specific competition. This 
equilibrium dynamics is also possible when the intensity of intra-specific competition intensity is less 
than that of inter-specific competition but under conditions that the maximum intrinsic growth rate of 
one species is not too large. The coexistence of two competing species occurs only if the system has 
four interior equilibria. We provide an approximation to the basins of the boundary attractors (i.e., 
the extinction of one or both species) where our results suggests that contest species are more prone to 
extinction than scramble ones are at low densities. In addition, in comparison to the dynamics of two 
species scramble competition models subject to Allee effects, our study suggests that (i) Both contest 
and scramble competition models can have only three boundary attractors without the coexistence equi¬ 
libria, or four attractors among which only one is the persistent attractor, whereas scramble competition 
models may have the extinction of both species as its only attractor under certain conditions, i.e., the 
essential extinction of two species due to strong Allee effects; (ii) Scramble competition models like 
Ricker type models can have much more complicated dynamical structure of interior attractors than 
contest ones like Beverton-Holt type models have; and (iii) Scramble competition models like Ricker 
type competition models may be more likely to promote the coexistence of two species at low and high 
densities under certain conditions: At low densities, weak Allee effects decrease the fitness of resident 
species so that the other species is able to invade at its low densities; While at high densities, scramble 
competition can bring the current high population density to a lower population density but is above the 
Allee threshold in the next season, which may rescue a species that has essential extinction caused by 
strong Allee effects. Our results may have potential to be useful for conservation biology: For example, 
if one endangered species is facing essential extinction due to strong Allee effects, then we may rescue 
this species by bringing another competing species subject to scramble competition and Allee effects 
under certain conditions. 

Key words: Allee effects, Scramble competition models, Contest competition models, Basins of 
attractions, Extinction, Coexistence 


1. Introduction 

An Allee effect is a biological phenomenon characterized by a positive correlation between a popula¬ 
tion density and its per capita growth rate at small population densities (Allee et al. 1949). A distinction 
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is made between a strong Allee effect and a weak Allee effect: a strong Allee effect refers as to a popula¬ 
tion that exhibits a “critical size or density”, below which population declines to extinction, and above 
which it may increase; while a weak Allee effect refers as to a population that lacks a “critical density”, 
but where, at lower densities, the population growth rate rises with increasing densities (Stephens et al. 
1999; Lidicker 2010). Allee effects have been detected in natural populations for a wide array of taxa 
and are believed to be strong regulators of extinction, influencing colonization success of invasive species, 
disease dynamics, and the long-run population demise in conservation. Mathematical models can help 
us understand these combined effects on persistence, expansion or extinction of biological species. There 
is a considerable amount of populations dynamics in the presence of Allee effects (e.g., Dennis 1989 & 
2002; Selgrade and Namkoong 1992; McCarthy 1997; Shigesada and Kawasaki 1997; Greene and Stamps 
2001; Keitt et al 2001; Fagan et al 2002; Wang et al 2002; Liebhold and Bascompte 2003; Schreiber 2003; 
Drake 2004; Zhou et al 2004; Petrovskii et al 2005; Taylor and Hastings 2005; Jang 2006; Aguirre et al 
2009; Egami 2009&2010; Thieme et al. 2009; Elaydi and Sacker 2010; Wang et al 2010; Liu et al 2011; 
Kang and Castillo-Chavez 2012; Kang and Castillo-Cliavez 2014a&b; Cai et al 2014; Peng and Kang 
2015; Livadiotis et al 2015) as well as various models in patchy environments (e.g., Amarasekare 1998a 
&1998b; Gyllenberg et al 1999; Ackleh et al 2007; Kang and Lanchier 2011; Kang and Armbruster 
2011; Kang and Castillo-Chavez 2014c). In another direction, evolutionary models with Allee effects 
have gained some attentions (Cushing and Hudson 2012; Kang and Udiani 2014; Cushing 2015). For 
example, Cushing and Hudson (2012) investigated the global dynamics of an evolutionary model for a 
population subject to strong Allee effects. One of their results suggests that evolution is beneficial in 
the sense that reduces the possibility of extinction due to an Allee effect. 

Resource competition can be defined in the contest-scramble spectrum competition according to re¬ 
source partitioning among competitors (Nicholson 1954). Contest and scramble-type competition are two 
extreme forms of competition, which are characterized by resource monopolization and resource sharing, 
respectively. Contest competition results in a constant number of survivors with enough resource gain 
against an initial density of competitors, while scramble competition (also called over-compensation) 
results in an increased number of survivors with decreased resource intake and an increased number of 
competitors (Calow et al. 1998).A well-known example for contest competition is the Beverton-Holt 
model while the well-known example for scramble competition is the Ricker model. Even though Allee 
effects are expected commonly in nature, they are often ignored in studies of species coexistence mech¬ 
anisms. Many species can experience different forms of competition at high densities and Allee effects 
at low densities with the consequence that their populations do not grow optimally at low densities 
and individuals compete with one another at high densities (Begon et al. 1996; Etiemme et al. 2002; 
Kang 2013). The potential effects of combining positive density dependences from Allee effects with 
negative density dependences coming from different types of competition are expected to generate rich 
dynamics due to the fact that outcomes of varied patterns of competition have different consequences for 
population dynamics and the evolution of individual traits (Godfray 1987; Lomnicki 1988; Ives 1989). 
For instance, with a species at a low density, in the invader state, it faces both intra-specific and inter¬ 
specific competition, but obtains a lesser effect of reduced intra-specific competition due to Allee effects 
compared with the standard investigations of models of competing species. As comparisons between 
intra-specific and inter-specific density dependences are critical to competitive coexistence, ignoring the 
potential that Allee effects greatly modify intra-specific competition is a critical oversight. For instance, 
Chesson and Ellner (1989) has pointed out the potential for major effects such as stochastic models 
may have interior invariant regions if low-density growth rates are depressed by Allee effects. Hopf et 
al (1993) also point out potential major implications for community structure of Allee effects uncovered 
by introducing them into a competition model. 

Recently, Kang and Yakubu (2011) and Kang (2013) have explored the population dynamics of two 
competing species when both species experience scramble competition (i.e., Ricker-type models) and 
Allee effects induced by predation saturations. The positive density dependence from weak Allee effects 
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can decrease the fitness of resident species so that the other species is able to invade at its low densities. 
Thus weak Alice effects can promote the permanence of two competing species at their low densities 
(Kang and Yakubu 2011). In the presence of strong Allee effects, scramble inter-specific competition 
can bring the current hight population density to a lower population density but being above the Allee 
threshold. Thus strong Allee effects can save species from potential essential extinction (Schreiber 2003) 
by bringing in scramble inter-specific competition (Kang 2013), therefore promote the existence of two 
competing species at their high densities under certain conditions. Understanding how different forms 
of competition in the presence of Allee effects affect population dynamics when species suffer from Allee 
effects, can advance our understanding of the extinction and establishment of species in ecological com¬ 
munities, with implications for conservation programs (Zhou et al 2004; Courchamp et al. 2009; Kang 
and Yakubu 2011; Kang 2013; Kang et al 2014 a&b). 

A recent work by Livadiotis and Elaydi (2012) considered a general framework of population models 
subject to strong Allee effects. We adopt and modify their definitions on mathematical formulation 
of strong Allee effects in single species models. Our main purpose is to investigate the dynamical 
outcomes in two species contest competition models when both species have Allee effects. The contest 
competition models in this article refer as to a generalized Beverton-Holt competition model while 
scramble competition models refer as to species that have one-hump growth functions (e.g., Ricker’s 
map) in its single state. More precisely, we would like to explore possible answers to the following 
ecological questions: 

1. What are the dynamical outcomes of a two species competition model when both species experience 
both strong Allee effects and contest competition? 

2. How do the dynamics of models with a strong Allee effect and contest competition compare to 
those of models with a strong Allen effect and scramble competition? 

3. What are the generic dynamical features of population models subject to strong Allee effects and 
different forms of competition at low and high densities? 


The paper is organized as follows: In Section 2, we define how we will mathematically formulate Allee 
effects, contest and scramble competition for single species population models and derive a generalized 
Beverton-Holt competition model that each species is subject to strong Allee effects. In addition, we 
provide the important results of single species models with strong Allee effects and contest competition 
that can be applied to two species competition models. In Section 3, we apply the theory of monotone 
dynamics to obtain sufficient conditions when our two species competition model has equilibrium dy¬ 
namics. This result can establish some generic dynamical features of contest competition models subject 
to Allee effects. We also give approximated basins of attractions of boundary attractors (i.e., initial 
conditions that lead to the extinction of one or both species). Moreover, we explore conditions that lead 
to the coexistence of two competing species. These results can help us understand the coexistence and 
the extinction conditions for two competing species that are subject to contest competition and Allee 
effects. In Section 4, we apply the theoretical results derived in Section 3 to a specific two competing 
species model with strong Allee effects and contest competition. Numerical simulations are performed 
to obtain more insights into the generic features of two species models with strong Allee effects and 
contest competition. The results obtained from Section 3 &4 can partially answer three questions pro¬ 
posed in the introduction. In Section 5, we discuss the dynamic outcomes of two competing species 
population models with different types of competition (i.e., contest and scramble competition) where 
each species has strong Allee effects. In addition, we summarize our results and mention some potential 
future studies. In Appendix, we provide detailed proofs of our theoretical results. 
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2. Model Derivations 


In this article, we are interested in the population dynamics of a generalized Beverton-Holt two 
species competition model subject to Allee effects. More specifically, we would like to explore how the 
interactions of contest competition and Allee effects affect the establishment and extinction of species. 
By understanding population dynamics of such models, we are able to obtain insights on the impact of 
different forms of competition combined with Allee effects on the persistence and extinction of species by 
comparing the dynamics to the case when two competing species both suffer from scramble competition 
and Allee effects. 

2.1. Single species population models 

First, we would like to introduce the terms of contest competition, scramble competition and Allee 
effects for one species population models rigorously. Let it t be a species population density in season f, 
then its population density in season t + 1 can be represented as the following equation: 

u t+ i = H(u t ) =u t h(u t ) (1) 

where H(u) > 0, u £ R+ is the growth function of species u and h{u ) is the per capital growth rate that 
is subject to the condition lim,^ h(u) = 0 which guarantees the boundedness of population (Kang and 
Chesson 2010; Kang 2013). We say that species u has contest intra-specific competition if lim,,^^ H(u) = 
c > 0 which implies that a constant number of survivors with enough resource gain against an initial 
density of competitors. We say that species u has scramble intra-specific competition if lim^^oo H(u) = 0 
which implies that an increased number of survivors with decreased resource intake and an increased 
number of competitors. Contest and scramble intra-specific competition are two extreme forms of 
competition that can result in very different consequences for population dynamics. For example, an 
important ecological example of contest intra-specific competition is the Beverton-Holt model (see Figure 
1(a)) whose H takes the form of H(u) = 1 ^ u while a well-known ecological example of scramble intra¬ 
specific competition is the Ricker’s model (see Figure 1(b)) whose H takes the form of H{u ) = ue r ~ au . 
The Beverton-Holt model has relatively simple population dynamics where population goes to 0 if r < 1 
and population goes to r — 1 if r > 1. While the Ricker’s model has very complicate population dynamics 
(e.g., periodic orbits, chaos) depending on the values of r. See Figure 1 for examples of single species 
models with contest intra-specific competition versus scramble intra-specific competition. 

For convenience, we define the following four conditions with respect to Model (1): 

Al: There exists some e > 0 such that > 0 for u £ (0, e). 

A2: h(u) > 0 and linq^oo h(u ) = a < 1. 

A3: There exists unique two numbers 0 < A < I\ , such that h(A) = h{K) = 1 and h'(A) > 1. In 

addition, h{ 0) < 1. 

A4: There exists a unique fixed point K , i.e., h{K) = 1. In addition, hfiS) > 1. 

Condition Al indicates that the per capita growth rate has a positive correlation with population density 
for small populations. Condition A2 indicates that the population is positive invariant (i.e., all future 
population is nonnegative for any nonegative initial population) and bounded. Condition A3 indicates 
that the population has two steady states A and K where A, called Allee threshold , is always unstable 
and K, called carry capacity , can be locally stable under certain conditions. Condition A4 indicates 
that species u is persistent and may be locally stable at its carry capacity K. Now we can define the 
following terms: 

• Species u has Allee effects if its population model (1) satisfies Condition Al, A2. 
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H(u): Contest Competition Model | |" 1 H(u): Scamble Competition Model 




(a) 


(b) 


Figure 1: Examples of single species models with contest competition (a) versus scramble competition (b). 


• Species u has strong Allee effects if its population model (1) satisfies Condition Al, A2, A3 (see 
blue curves in Figure 2-3). 

• Species u has weak Allee effects if its population model (1) satisfies Condition Al, A2, A4 (see 
black curves in Figure 2-3). 

Our definition of Allee effects is more general than Livadiotis and Elaydi (2012) since they only define 
strong Allee effects as the following three conditions: i) h'(u) > 0 for u £ (0,e) for some e > 0; ii) 
h( 0) < 1 and iii) There exists a unique K > 0, such that h{K) = 1, h'(K) < 0. These two definitions can 
be considered the same. When species u is subject to strong Allee effects , i.e., its population model (1) 
satisfies Condition Al, A2, A3, its population dynamics can be summarized as follows: When initial 
population density is below A, the population goes extinct; When initial population density is above 
A , the population may be able to sustain; when h( 0) is large enough, then (1) may have the essential 
extinction (i.e., for any positive initial condition, the population converges to 0 with probability 1) due 
to strong Allee effects (Schreiber 2003). See Figure 2-3 for an illustration of strong Allee effects (blue 
curves in both figures) versus weak Allee effects (black curves in both figures). 

Outcomes of varied patterns of intra-specific competition (e.g., contest competitions or scramble 
competitions) have different consequences for population dynamics (1) when species experiences Allee 
effects. A well-known contest competition model introduced by Thomson (1993) for fisheries is called 
the Sigmoid Beverton-Holt model that can be described as 

ru S ~ 1 

u t +i = H(u t ) = u t h(u t ) = u t -—-—t (2) 

1 + Ut 

where r,S > 0. Model (2) is a depensatory generalization of the Beverton-Holt stock-recruitment rela¬ 
tionship used to develop a set of constraints designed to safeguard against overfishing (Thomson 1993; 
Myers et al 1995; Myers and Mertz 1998; Stoner and Ray-Culp 2000; Myers 2001; Gascoigne and Lipcius 
2004; Harry et al 2012). One important feature of (2) is that it exhibits the Allee effects if S > 1. The 
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~~ H=u No Allee effects ~~ Weak Allee effects ~~ Strong Allee effects 



Figure 2: An illustration of strong Allee effects versus weak Allee effects of the growth function H(u) = uh(u)\ The red 
curve is H = u or h = 1; the green curve is the case when (1) has no Allee effects; the blue curve is the case when (1) has 
strong Allee effects and the black curve is the case when (1) has weak Allee effects. 


dynamical properties of (2) have been proposed by Harry et al (2012). Here we summarize the important 
dynamics of (2) that will be used later for our two-species model as the following proposition: 

Proposition 2.1 (Dynamical of (2)). Let r, 5 > 0 and (2) has an initial condition Xq > 0. If S > 1, we 
define r cr a = 5(5 — l)" 1 . Then the following statements are true: 

1. If 5 £ (0,1), then Model (2) has two non-negative equilibria: an unstable equilibrium 0 and a stable 
positive equilibrium I\. 

2. If 5 > 1 and r < r cr it, then the only equilibrium of Model (2) is 0 and it is locally stable. 

3. If 5 > 1 and r > r cr it, then Model (2) has three non-negative equilibria: 0, A and K such that 
0 < A < r(<5 ~ 1} < K. The 0 equilibrium is stable with the basin of attraction [0, A); A is unstable 
(repellor), while I\ is stable with the basin of attraction (A, oo). 


Notes: Harry et al (2012) gave similar results of Proposition 2.1 as their Proposition lbut without 
rigorous proof. Here, we provide a detailed proof in the Appendix. Proposition 2.1 implies that Model 
(2) has Allee effects for any r, <5 > 0. Precisely, Model (2) has weak Allee effects if 5 £ (0,1) while Model 
(2) has strong Allee effects if 5 > 1 and r > r cr it■ In addition, we would like to point out that when 
5 = 2 , Elaydi and Sacker (2009) have extended (2) to the following generalized form: 


u t +i 


u t 


du t + e 
u\ + bu t + c 


( 3 ) 


where b is called the shock recovery parameter and c can be considered as a carrying capacity. The 
dynamics of (3) is similar to (2) when 5 = 2. Define a single species population model as follows 


Ut+l 


H(u t ) = u t h(u t ) = u t 


ru 


<5-1 


a + Uf 


( 4 ) 
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- h=l 

No Allee effects 

Weak Allee effects — 

— Strong Allee effects 



Figure 3: An illustration of strong Allee effects v.s. weak Allee effects of the per capita growth rate h(u) = 
:corresponding to its growth function H(u) in Figure 2: The red curve is H = u or h = 1; the green curve is 
the case when (1) has no Allee effects; the blue curve is the case when (1) has strong Allee effects and the black curve is 
the case when (1) has weak Allee effects. 


where r is the maximum intrinsic growth rate and a can be considered as a carrying capacity as the 
parameter c in Model (3). Then the dynamics of Model (4) can be obtained as the following corollary 
by directly applying the results from Proposition 2.1 if we define r crit = a s S(S — l)* -1 : 

Corollary 2.1 (Dynamical of (4)). Let r, a, 5 > 0 and (4) has an initial condition Xq > 0. If S > 1, 
we define r" rit = a 1 ^ s 6(S~ l)? -1 . Then the results of Proposition 2.1 can directly apply to Model (4). 
In particular, if a, 6 > 1 and r > rf rit , then Model (4) has two interior equilibria A a and I\ a such that 
0 < A < A a < K a < K where A , K are defined in Proposition 2.1. 

The scramble intra-specific competition models subject to Allee effects can take a form of the Ricker’s 
model, for example, 


u t +i = u t h(u t ) = u t e r(1 Ut) 1+i "f (5) 

where r represent the intrinsic growth rate (also the scramble intra-specific competition coefficient); m 
represents predation intensities and b represents the product of the proportional to the handling time 
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and the carrying capacity. For Model (5), we have 


h(u) = e r( ' 1 1 + fc ' u and h'(u) 


[6w-r(l + 6u) 2 ]e r(1 - M) ^Tfk 
(1 + bu) 2 


Thus, according to the definitions above, Model (5) has Alice effects if mb > 0; Model (5) has weak Allee 
effects if bm > r > m; and Model (5) has strong Allee effects if r < m < ri ' 1 ^ , b > 1. In addition, 
Model (5) has more complicate dynamics (e.g., chaos, essential extinctions) than the contest competi¬ 
tion model subject to Allee effects (2). The detailed dynamics of (5) has been studied by Schreiber (2003). 


Another example of scramble competition with Allee effects is a generalized Maynard Smith-Slatkin 
model that can be described as follows: 

ruf- 1 

ut+ i = H(u t )=u t h(u t ) = u t —— l d (6) 

1 + u° + buf 

where all parameters are positive; r, 8 have the same biological meanings as in Model (2) and b , d have 
the same biological meaning as in Model (3). Model (6) can exhibit similar complicated dynamics as 
Model (5) due to the generic feature of scramble competition. Here we summarize the basic dynamics of 
(6) that will be used in the later sections as the following theorem: 

Theorem 2.1 (Dynamics of Model (6)). If 8 > 1, we define r cr u as follows 

( 8 y i_,5)/d \ d t ( 8 \ S/d ' 

Tcrit ~ \b(d-8)) d — 8 + \b(d - 8)) 

Depending on the values of r, b , < 5 , d, the basic dynamics of (6) can be classified as the following cases: 

1. If 8 > d, then the trajectory of Model (6) with initial condition Xq > 0 converges to one of its 
equilibria. More specifically, we have the following two cases: 

(a) If 8 < 1, then Model (6) has two non-negative equilibria 0 and K where 0 is unstable and K 
is globally stable. 

(b) If 8 > 1, then u = 0 is unstable and Model (6) can have 0 or two or even four positive equilib¬ 
ria. Moreover, Model (6) is globally stable at 0 when it has no positive interior equilibrium. 

2. If 8 = d, then Model (6) is reduced to Model (2) where its dynamics has been stated in Proposition 
2 . 1 . 

3. If 8 < d, then u = 0 is always locally asymptotically stable. In addition, 

(a) If 8 < 1, then Model (6) has two non-negative equilibria 0 and K where 0 is unstable. 

(b) If 8 > 1 and r > r cr u, then Model (6) has at least two distinct positive roots A, K such that 
0 < A < u c < K. 

f. If both 8 and d are positive integers, then Model (6) has at most two distinct positive interior 
equilibria. 

Notes: Theorem 2.1 implies that Model (6) can exhibit strong Allee effects if <5 > maxjl, d} or 1 < 8 < 
d,r > r cr it- Model (2) can exhibit strong Allee effects only if 8 > 1 and r > r cr it- Even though Part 1 (b) 
and Part 3 (b) of Theorem 2.1 do not state that Model (6) has exactly two interior equilibria, numerical 
simulations suggest that Model (6) has at most two interior equilibria. In addition, if 1 < 8 < d and 
r > rcrit ? then Model (6) is a unimodal map that can have very complicated dynamics as the Ricker’s 
map including the essential extinction due to strong Allee effects. Numerical simulations suggest that 
increasing values of r, b, 8 can destabilize the system while increasing value of d can stabilize the system 
(see Figure 5-4 as examples). 











Figure 4: Periodic time series of Model (6) when 5 = 2; d = 5; b = 0.1 and r = 7.5. 


2.2. A generalized Beverton-Holt competition model with Alice effects 

A generalized Beverton-Holt two species competition model with Allee effects from single species 
models (2) and (6) can be represented by the following two equations: 


x t +i 


Vt+i 


rixj 1 

l+x 5 f+ b lV s t 3 

yf 2 

1 + y S t *+b 2 X S t * 


( 7 ) 

( 8 ) 


where rj ,i = 1,2 describes the maximum intrinsic growth rate of each species; b^i = 1,2 describes 
the relative inter-specific competition coefficient to its intra-specific competition; = 1,2 describes 
the nonlinearity of the intra-specific competition, while 8i,i = 3,4 describes the nonlinearity of the 
inter-specific competition. Define Condition HI: 

HI: Si > 1, Ti > r' crit , i = 1,2 where r l crit = Si(Si - 1)A -1 ,* = 1,2. 

According to Proposition 2.1, each species of Model (7)-(8) suffers both contest intra-specific competition 
and strong Allee effects if Condition HI is satisfied. Due to our interests, we study the population 
dynamics of (7)-(8) when Condition HI is satisfied for most cases. 


3. Dynamical properties of a generalized Beverton-Holt competition model with Allee 
effects 

In this section, we explore the basic dynamical property of Model (7)-(8) in terms of its monotonicity, 
the extinction and coexistence of two competing species. If Condition HI is satisfied, Model (7)-(8) 
always has the following five boundary equilibria: 

Eg = (0,0), E Xa o = (Ai,0), E Xk o = (Ai,0), Eo Va = ( 0 ,^ 2 ) and Eo yk = (0,/\2) 
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Chaotic dynamics 



Figure 5: Chaotic time series of Model (6) when 5 = 2;d = 5;6 = 0.1 and r = 25. 


where Ai,i = 1,2 is the Allee threshold of species i and Ki,i = 1,2 is the carrying capacity of species i, 
thus, 0 < Ai < K i: i = 1,2. Simple calculations indicate that E 0 ,E Xk0 and E 0yk are locally asymptot¬ 
ically stable, thus these three boundary equilibria are also attractors for Model (7)-(8). This gives the 
following lemma: 

Lemma 3.1 (Basic properties). With respect the map defined by Model (7)-(8) the positive cone R+ is 
forward invariant and the rectangular [0, r i] x [0, r<f\ is absorbing. If Condition HI is satisfied, then it has 
five boundary equilibria Eq, E Xa o, E Xk o, Eoy a and Eo Vk where Eq, E Xk o and Eo yk are locally asymptotically 
stable and E Xa0 ,E Va0 are unstable. In addition, let Xi = {(cc, 0) £ R)j_} and X 2 = {(0,y) £ R(j_}, then 
both Xi, i = 1,2 are positive invariant sets for Model (7)-(8). 

Notes: It is easy to check that the set [0, rf\ x [0, r 2 ] is forward invariant and all orbits in the positive 
cone R^ lie in this set after one time step. The stability of the boundary can be obtained by straight 
forward calculations, i.e., calculating the eigenvalues of the Jacobian matrix evaluated at the boundary 
equilibria. Thus, we omit the details here. Also, we would like to point out that E Xa o and Eo Va are 
saddle nodes: E Xa o is a saddle node with the unstable manifold lying on the x-axis and the stable 
manifold tangent to the vertical eigenvector emanating from E Xa0 ; similarly, E 0ya is a saddle node with 
the unstable manifold lying on the y-axis and the stable manifold tangent to the horizontal eigenvector 
emanating from Eo Va . 

3.1. Competitive systems 

Monotonicity properties have been used successfully to study the dynamical behavior of solutions 
to differential equations (e.g., see [Smith 1995] for an overview) and difference equations (e.g., Selgrade 
and Ziehe 1987; Dancer and Hess 1991; Hess and Lazer 1991; Smith 1998; Wang and Jiang 2001; Clark 
et al 2003). A competitive map is one type of monotone systems. 

To continue our study, we first define a partial order on R 2 so that the positive cone in this new par¬ 
tial order is the fourth quadrant. Specifically, for u = (ui,u 2 ),v = (vi,v 2 ) £ R+, we say that u <k v if 


10 





u\ < Vi and i>2 < «2- Two points u, v £ R+ are said to be related (or iv-ordered) if u <k v or v <k u. A 
strict inequality between points can be defined as u <k v if u <k v and u ^ v. A stronger inequality can 
be defined as u «k v if u\ < Vx and V 2 < U 2 - A map T : IntM. 2 + —> Int is competitive (or AT-order 
preserving) if u <k v implies that T[u) <k T(v) for all u,v £ Int R+. A map T : IntM 2 , —> IntM 2 + is 
strictly competitive (strongly competitive) if u <k v implies that T(u) <k T(v) (' T(u ) «k T{v)) for all 
u,v £ IntM}_. Clearly, being related is an invariant under iteration of a strongly competitive map. Let 
J be the Jacobian matrix of a map T, then J is A'-positive (A'-strongly positive) if its diagonal entries 
are nonnegative (positive) and its off-diagonal entries are nonpositive (negative). 


Theorem 3.1 (Competitive systems case 1). If S 1 S 2 > 6364 , then Model (7)-(8) is strongly compet¬ 
itive and all orbits in the positive cone converge to an equilibrium. 

Notes: The proof of Theorem 3.1 use the monotone properties of Model (7)-(8), which is provided in 
the Appendix. This theorem indicates that Model (7)-(8) has equilibrium dynamics (e.g., no nontrivial 
periodic orbits) whenever A > ^ 3 ^ 4 - Recall that Si, i = 1,2 measures the intensity of the intra-specific 
competition of species, and Si,i = 3,4 represents the inter-specific competition. Then we can define 
as intensities of the inter- and intra-specific competition, respectively. Then the result is that 
the dynamics equilibration when the intensity of intra-specific competition intensity exceeds that of 
inter-specific competition. 

From the proof of Theorem 3.1, we can see that Model (7)-(8) is strictly competitive in and is 
strongly competitive in Inf R^_. For the case that 6362 < J 3 J 4 , we have the following theorem: 


Theorem 3.2 (Competitive systems case 2). Assume that 6362 < S3S4. If 

( 2 \ 1/Si ( J1J2 \ 1/S3 

^62^3^4 — <^1^2]/ 2 \6 i [ c )3&1 — S182] J 

then every orbit of Model (7)-(8) converges to an equilibrium. 

Notes: The proof of Theorem 3.2 use the monotone properties of Model (7)-(8) as well as the prop¬ 
erties of non-invertible maps. The detailed of proof has been given in the Appendix.In the case that 
J 1 J 2 < J 3 J 4 , Theorem 3.2 indicates that the dynamics of (7)-(8) can still have equilibrium dynamics 
as the case when 8162 > J 3 J 4 (Theorem 3.1) if the maximum intrinsic growth rate of one species is 
not too large. However, when this additional condition is not satisfied, Model (7)-(8) can have non¬ 
equilibrium attractors as suggested by results of Terescak (1996) and Smith (1998). For example, if we 
let <5i = S 2 = 8,63 = 84 = d,r 1 = r 2 ,bi = 62 , then the dynamics of (7)-( 8 ) on the invariant manifold 
H x =y = {( x,y ) £ R+ : x = y} can be reduced to Model ( 6 ) which can have complicated dynamics. See 
Figure 4-5 as examples for different dynamics on fl x=y when 8 = 2;d = 5; 6 = 0.1 and r = 7.5 (Figure 
4) and r = 25 (Figure 5). 


3.2. Basins of attractions of extinction equilibria 

Assume that Condition HI is satisfied for Model (7)-(8). 

O 0 = [ 0 ,Ai] x [0,A 2 ], O ex = {{x,y) £ R+ : x < Ax} 

and Condition H2-H3 as follows: 

H2: $1 > 1 and r\ > rff it where r a c f it = {ai) 1 / Sl 8 x{ 8 i - 

H3: 82 > 1 and r 2 > rff it where r“^ t = (c^) 1 ^^^ — 


We can define the following sets: 
and O ey = {(a;, y) £ R+ : y < A 2 } 

l )^ -1 and ax = 1 + bx{A 2 ) S3 . 
l )^ -1 and 02 = 1 + b 2 {Ax) 5i . 
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Condition H2 or H3 indicates that species suffers from strong Allee effects. Assume that Model (4) 
satisfies Condition H2, then according to Corollary 2.1, it has two interior equilibria A® 1 and K where 
AI 1 is the Allee threshold and Kf 1 is the carrying capacity. Similarly, if Model (4) satisfies Condition 
H3, then it has two interior equilibria A® 2 and Kf 2 where A^ 2 is the Allee threshold and K% 2 is the 
carrying capacity. Let 

O x = {(x,y) £ R+ : x > A® 1 , y < A 2 } and O y = {(x,y) £ M+ : x < A 2 , y > Af 2 }, 
then we can show the following theorem: 


Theorem 3.3 (Basins of attractions of boundary attractors I). Assume that Model (7)-(8) satisfies 
Condition HI , then 

• For any initial value in O ex , we have Hindoo xt = 0. 

• For any initial value in O ey , we have lim^oo y t = 0. 

• For any initial value in Oq, we have lim t _ ) . 00 (cct, yt) = Eq. 

If Model (4) satisfies Condition H2, then for any initial value in O x , we have 

lim {x t ,y t ) = E x o- 

t—¥ OO 1 

If Model (4) satisfies Condition H3, then for any initial value in O y , we have 

lim {x t ,yt) = E 0y 

t—> OO 2 

Notes: Theorem 3.3 (see the proof in the Appendix) provides an approximation of the boundary at¬ 
tractor’s basins of attractions for Model (7)-(8) when each species suffers from contest competition and 
strong Allee effects. In the absence of other species, each species’ extinction region is [0,Aj] while in 
the presence of other species, the extinction region becomes larger due to the inter-specific competition, 
thus, competition cannot save species from extinction at their low densities. In order to maintain species 
x’s population, its initial condition should be large enough and the other species population should be 
small enough, e.g., the initial condition should be larger than A ® 1 > A\ while species y’s initial condition 
should be less than A 2 - In addition, Theorem 3.3 implies that the basins of attractions of the boundary 
attractors E Xki0 and E 0yk2 are unbounded (see Figure 12, 14 and 16). 

Define 

° l ex = Ufc> 0 {(*. kx ) e IntB.% : X > (^l 1 ) ' 3_Sl } 

° l ey = U>0 {{ k y,y)£ I nm\:y> 

and let O l 0 = O l ex n O l ey . Then we have the following theorem: 

Theorem 3.4 (Basins of attractions of boundary attractors II). Assume that Model (7)-( 8 ) satisfies 
Condition HI. 

• If 63 > Si, then for any initial value in O ex U O l ex , we have lim^oo x t = 0. 

• If 64 > 5 2 , then for any initial value in O ey U O l ey , we have lim^oo y t = 0. 

• If S 3 > 8 1 and 84 > S 2 , then for any initial value in Oq U Oq, we have lim t _ > . 00 (a;t, yt) = Eq. 
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Notes: The detailed proof of Theorem 3.4 is provided in the Appendix. This theorem indicates that if a 
species’ relative competition degree is less than 1 (i.e., S 3 > <5i or £4 > <$ 2 ), then its extinction region (i.e., 
the initial conditions that lead to the extinction of two species) consists of two distinct components, of 
which one is stated in Theorem 3.3 and the other one is the area that when its competitor’s population 
is large enough. This phenomenon is is a typical property of non-invertible maps which is caused by the 
large degree of the inter-specific competition, i.e., £3 > <5i or £4 > <5 2 . Theorem 3.4 also indicates that 
the basins of attractions of Eq may be unbounded. In fact, the basins of attractions of Eq is a connected 
and unbounded region when there is no interior equilibrium (see Figure 11 and 12). While the basins 
of attractions of Eq consists of two parts when there is an interior equilibrium: one part is a bounded 
region including the neighborhood of E 0 while the other part is a connected unbounded region including 
very large initial values of both species (see Figure 14 and 16). 


3.3. Persistent equilibria 

Assume that Model (7)-(8) satisfies Condition HI, then according to Theorem 3.3, it always has 
three boundary attractors. This suggests that the coexistence of two species is possible only if there 
exists an interior equilibrium. We are going to prove the following theorem regarding the number of 
interior equilibria. 

Theorem 3.5 (Interior equilibria). Assume that Model (7)-(8) satisfies Condition HI. Let 


Fi(x) 



i/S 3 

and A 2 (j/) 



1/54 


then we have the following statements: 

• If Fi(x c ) < A 2 or F 2 (y c ) < A\, then Model (7)-(8) has no interior equilibrium. 

• If A 2 < Fi(x c ) < K 2 , K\ < F2(y c ) or A\ < F 2 (y c ) < Ad, K 2 < F\(x c ), then Model (7)-(8) has 
two interior equilibria. 

• If Fi(x c ) > A ' 2 and A 2 (y c ) > K\, then Model (7)-(8) has four interior equilibria, 
where x c = ri ^ ^ and y c = r2 ^ 2 1 ' 1 . 

Notes: See Figure 6-10 regarding the generic nullclines of Model (7)-(8). Theorem 3.5 (its proof stated 
in the Appendix) gives us two scenarios when Model (7)-(8) has interior equilibria where coexistence 
occurs only if Model (7)-(8) has four interior equilibria (see Figure 15 and Figure 16). Sufficient conditions 
guaranteeing four interior equilibria is that Fi{x c ) > K 2 and A 2 (y c ) > Ad, i.e. 


x c < 


nx. 


Si 


1 + xt 1 + b\Kl 


and y c < 


X 2 yj 2 

1 + y s c 2 + b 2 K[ i 


which indicates that any initial condition in [x e ,A'i] x [y c ,AT 2 ] lead to a locally asymptotically stable 
interior equilibrium according to the monotonicity of the system. Biologically, this implies that both 
species are able to persist under proper initial conditions if Fi(x c ) > A " 2 and A 2 (y c ) > Afi holds. In the 
next subsection, we explore the local stability of a symmetric system. 


S.f. A symmetric competition model with Allee effects 

In order to obtain more insights on the coexistence of two competing species, we focus on a symmetric 
case of (7)-(8) where 7*1 = r 2 = r, 61 = b 2 = b , = S 2 = 5, S 3 = d 4 = d which gives the following model 


rx t 

1 + xf + byf 

(9) 

ry s t 

1 + yt+bxf 

(10) 
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No Interior Equilibrium 



Figure 6: Schematic nullclines when Model (7)-(8) has no interior equilibrium of Case 1, i.e., 
curve is F\(x) = ^ rix 1 — b ~ x 1 -1 ^ where its maximum Fi(x c ) occurs at x c = - 1 ^ ; 

F2(y) = ( r2yS2 ~ 1 b - yS2 


1 /Hi 


where its maximum F -2 (y c ) occurs at y c = 7-2 A 2 —i). 


F 2 (i ic) < Ai: The red 
While the blue curve is 


Understanding the dynamics of (9)-(10) also can help us obtain a better understanding of the dynamics 
of asymmetric cases due to the robust property of our system (7)-(8). 

It is easy to verify that (9)-(10) is reduced to the 1-D model described by (6) if Xq = yo, i.e., the 
manifold defined by £l y = x = {(a;, z/) £ : x = y} is invariant. According to Theorem 2.1, Model 

(9)-(10) has at least two symmetric equilibria (A, A) and {K,K) such that 0 < A < < x if 


d > 5 > 1 and r > r cr i t 


6 

b(d - 6) 


(1 S)/d 


d 

d-S 


f s y /d ' 

{b(d-S)J 


( 11 ) 


Now we have the following theorem regarding the stability of (K, K). 

Theorem 3.6 (Stability of the symmetric equilibrium). Assume that Model (9)-(10) satisfies the 
inequalities (11). Then (9)-(10) has at least two symmetric equilibria (A, A) and (K,K) such that 0 < 

A < < x. Moreover, the symmetric interior equilibrium ( x*,x*) is locally asymptotically 

stable if 


5 + b(6 - d){x*) d S + b(S + d){x*) d 

1 + (x*) 5 + b(x*) d < an 1 + (x*) s + b(x*) d 

While the symmetric interior equilibrium ( x , x) is unstable if 


6 + b(S-d)(iT*) d S + b(S + d){x*) d 

1 + (x*) s + b(x*) d > ° r 1 + (x*) s + b(x*) d 


( 12 ) 


(13) 
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No Interior Equilibrium 



Figure 7: Schematic nullclines when Model (7)-(8) has no interior equilibrium of Case 2, i.e., F\(x c ) < A 2 : The red 
curve is Fi(x) = 


\ 1 /S 3 / r- , \ 

) where its maximum Fi(x c ) occurs at x c = ' 1 ^ While the blue curve is 


F 2 (y) = ( 


_ I r 2 y 2 -tr 2 -l 
62 


1 /Si 


—) 

where its maximum (y r ) occurs at y c = 7-2 A 2-11 . 


Notes: Theorem 3.6 indicates that it is impossible for the interior equilibrium of the symmetric model 
(9)-(10) to have Neimak-Sacker bifurcation from the invariant Q x=y since the eigenvalue that determines 
the points moving towards or away from Q x=y is always positive (see the detailed proof in the Appendix). 
This suggests that (9)-(10) may have relatively simple dynamics. 


4. Application to an example 

In this section, we let 5\ = 82 = 2, then Model (7)-(8) can be rewritten as follows: 


x t +i = 


yt+i = 


nx 2 t 


-5a 


1 + A + hy t 

riyt 

1 + vt + 


(14) 

(15) 


In the absence of species y, i.e., yo = 0, we have (14)-(15) is reduced to a single species model (2) when 
r = ri and 8 = 2, i.e., 


x t +i = 


ri xj 
1 + xf/ 


(16) 


Model (16) is a well-known fisheries model when fish is subject to contest intra-specific competition and 
strong Allee effects. The dynamics of (16) is very simple: 1. If ri < 2, then (16) converges to 0 for any 
Xq > 0. 2. If ri > 2, then (16) has two interior equilibria A\ and K\ where 

A, = 71 and Kl = 1 + FrTyl 
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Two Interior Equilibria 



Figure 8: Schematic nullclines when Model (7)-(8) has two interior equilibria of Case 1, i.e., K 2 < Fi(x c ), A\ < F 2 (y c ) < 


K\\ The red curve is F\(x) = 
curve is F 2 (y) = (^ r2V 2 — b ~ y 


) i/83 / f •< \ 

where its maximum F\ (x r .) occurs at x c = ri - j - —While the blue 


where its maximum (y r ) occurs at y c = ' 2( ^ —i). 


such that (16) converges to 0 if cco < A± while (16) converges to Ki if Xq > A 1 . 

Similarly, in the absence of species x, i.e., x 0 = 0, (14)-(15) is reduced to a single species model (2) 
for the case that r = r\ and S = 2, i.e., 


Vt+i 


F2 Vt 

1 + 2/t 


where it also has two interior equilibria A 2 and K 2 where 


(17) 


^-2 


r 2 ~ V r 2 - 4 

2 


and K 2 


t*2 + \]r\ - 4 
2 


Due to our interests, we focus on the case that r t > 2 ,i = 1,2 for Model (14)-(15). Then apply 
Theorem 3.1, Theorem 3.3, Theorem 3.2, Theorem 3.4 and Theorem 3.5, we obtain the following corollary 
regarding the dynamics of (14)-(15): 

Corollary 4.1 (Dynamics of (14)-(15)). Assume that Ti > 2. Define 


r 


a 1 

crit 




and r a c f it 




l/«4 


Then the following statements are true: 

1. If 6364 < 4, then every orbit of Model (14)-(15) with any initial condition in converges to one 
of its equilibria. 
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Two Interior Equilibria 



Figure 9: Schematic nullclines when Model (7)-(8) has two interior equilibria of Case 2, i.e., A 2 < Fi(x c ) < K 2 , K\ < 


FiiVc)'- The red curve is F\(x) = (’ 
- 1 _ 

&2 


\ ±/03 r (5 1) 

-J where its maximum F± (x c ) occurs at x c = ' 1 ^ While the 


/ ^2 — 1 ^2 1 \ 1 ^ ^ 

blue curve is i* 2 (y) = ( — --, ~ v — = — I where its maximum F -2 (y r ) occurs at y c = ’ 2 *■ ? ~ ' . 


2. If ( 53<54 > 4 and 


r i < 


1/(54 


or r 2 < 


4 \ 1/53 

bi [( 53^4 — 4 ]) 


— 4] j 

then every orbit of Model (14)-(15) with any initial condition in R^ converges to one of its equi¬ 
libria. 

3. Model (14)-(15) always has three boundary attractors Eq, E Xki o, Eo Vko where 
E„ = (0,0), = ( n + f^ ,0) and = (0, 

Moreover, 

• For any initial value in O ex , we have Hindoo x t = 0. While if 83 > 2, then for any initial 
value in O ex U O l ex , we have lim^oo Xt = 0. 

• For any initial value taken in O ey , we have lim^oo y t = 0. While if f 84 > 2, then for any 
initial value in O ey U O l ey , we have Hindoo y t = 0. 

• For any initial value taken in Oo, we have lim^oo (xt, yt) = Eq. While if 83 > 2 and 84 > 2, 
then for any initial value in OqKJO 1 0 , we have limt_>. 00 (xt, yt) = Eq. 

• If r\> rff it , then for any initial value in O x , we have 

lim (x t ,y t ) = E x 0 . 

t—¥ OO 1 

While if r 2 > then for any initial value in O y , we have 


t l_im (x t , y t ) = E 0yk2 


17 






















Four Interior Equilibria 





Figure 10: Schematic nullclines when Model (7)-(8) has four interior equilibria, i.e., F\{pc c ) > K 2 and Fztyc) > The 
red curve is F±(x) = ^ — 

F*(y) = ( rag<a -fc 2 - /a - 1 


&1 
1 /S* 


1/^3 ^ 

- J where its maximum F\(x c ) occurs at x c = ' 1 yWhile the blue curve is 


where its maximum Fhfec) occurs at y c = r2 ^~ 1 ^ . 


4 . Model (14)-(15) has no interior equilibrium if 


r \! 4- l \ 3 ^ r 2 - V ?, 2 -4 r l/ 4 — 1 

bi ) 


< 


l/<5 4 


< 


n 


- \/^- 4 


Model (14) -(15) has two interior equilibria if 

7*2 — x /^2 — 4 ( r l /^ - 4 \ 1 ^ 3 ^2 + \/r| - 4 


< 


61 


< 


and 


"f/4 - 1^ 1 ^ 4 ^ n + \Jr\ - 4 


ri - - 4 < /r|/4^\ 1/54 < rq + -4 ^ / r^/4- l \ 1/<53 > r 2 + \/r| - 4 

2 \ b 2 J 2 \ bi J 2 


Model (7)-(8) has four interior equilibria if 
,2/4 _ 1X^3 r2 + y/|34 


> 


and 


(^) 


1/54 


> 


ri + V^i ~ 4 


The corollary above does not give us information on the dynamic patterns of the system. Due 
to the complexity of the model, we use numerical simulations to investigate their dynamical patterns. 
According to simulations, we summarize its dynamical patterns as follows: 

• No interior equilibrium: When Model (14)-(15) has no interior equilibrium (see Figure 11), 
it has exactly three boundary attractors, i.e., (0,0), (a; 2 ,0) and (0, 2 / 2 ), and no interior attractor. 
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The basins of attractions of these three attractors are unbounded (shown in Figure 12) where the 
extinction area (i.e., basins of attractions of Eq = (0,0)) is an open and unbounded region in 
red. This implies that the inter-specific competition make both two competing species prone to 
extinction. 


Nullclines 



Figure 11: Nucllines of (14)-(15) when n = 7*2 = 2.5, 61 = 62 = 15 and < 53=54 = 5. 

• Two interior equilibrium: When Model (14)-(15) has two interior equilibrium (see Figure 13), 
it has exactly three boundary attractors, i.e., (0,0), (x 2 , 0 ) and ( 0 , 2 / 2 )) and no interior attractor. 
The basins of attractions of these three boundary attractors are unbounded (shown in Figure 14): 
The basins of attractions of (x2,0) and ( 0 , 2 / 2 ) are connected and unbounded while the basins of 
attractions of Eq consists of a bounded-connected and an unbounded-connected region. 

• Four interior equilibrium: This is the only case when Model (14)-(15) has an interior attractor 
(i.e., the coexistence of two species). When Model (14)-(15) has four interior equilibrium (see 
Figure 15), it has four attractors of which three are boundary attractors: (0,0), (X 2 , 0) and (0, 2 / 2 ), 
and one is an interior attractor (x*,y*) whose basins of attractions is a bounded and connected 
region. The basins of attractions of these four attractors are shown in Figure 16 where the features 
of basins of attractions of boundary attractors are the same as the two interior equilibria case. 

The observations above for Model (14)-(15) should be able to apply to a general model (7)-(8) due to the 
robustness of the model. In addition, we would like to point out that Model (7)-(8) can have coexistence 
of two species only if it has four interior equilibria. 
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Basins of Attractions 



Initial conditions of species x 


Figure 12: Basins of attractions of (14)-(15) when r\ = 7*2 = 2.5, 61 = 62 = 15 and £3 = £4 = 5: The black region is 
the basins of attractions of ( 0 ,?/ 2 ); The cyan region is the basins of attractions of (# 2 , 0 ); The blue region is the basins of 
attractions of the interior attractor which is an interior equilibrium and the red region is the basins of attractions of ( 0 , 0 ). 


5. Discussion 

Competition is an important ecological process which, in the short run, can cause a reduction in the 
number of species living within an area, by preventing very similar species from co-occurring; while in 
the long term, is likely to increase species diversity, by acting as a force for specialization and divergence. 
Competition can vary between two extreme forms: contest versus scramble. Allee effect is an another 
important ecological process that leads to a positive correlation of the per capita growth rate of a species 
and its low population density. The interplay of Allee effects and intra-specific competition can result 
in a critical threshold, called Allee threshold, below which species goes to extinction; while above which 
species may persist. This is referred to strong Allee effects. When there is no such threshold, it is referred 
as to weak Allee effects. Many species in nature are subject to both Allee effects and competition. The 
combinations of Allee effects and different type of competition can generate varied dynamical outcomes, 
especially for the coexistence of species. 

Including Allee effects in competition models can have huge potential impacts on diversity and com¬ 
munity structure (Chesson and Ellner 1989; Hopf et al 1993). Allee effects impose a cost of rarity 
(Mechanism R). When communities of organisms are subject to a cost of rarity imposed by Allee effects, 
they should be comprised of distinct species with differences in resource use that are greater than ex¬ 
pected by chance (Hopf et al 1993). Kang and Yakubu (2011) and Kang (2013) studied the population 
dynamics of a two competing species model where each species suffers from Allee effects and scramble 
competition: Kang and Yakubu (2011) showed that weak Allee effects may promote the permanence of 
two competing species at their low densities. The biological explanation for this is that weak Allee effects 


20 








































































































































































Figure 13: Nucllines of (14)-(15) when n = r 2 = 2.5, fei = &2 = -05 and <53 = <54 = 5. 

decrease the fitness of resident species so that the other species is able to invade at its low densities. 
Kang (2013) showed that strong Allee effects may save two competing species from essential extinctions 
at their high densities. The biological explanation for this is that scramble competition can bring the 
current high population density to a lower population density but is above the Allee threshold in the next 
season with the consequence that both competing species are able to persist. To extend the study by 
Kang and Yakubu (2011) and Kang (2013), this article studies the population dynamics of a generalized 
Beverton-Holt two species competition model where each species has strong Allee effects to explore the es¬ 
tablishment and the extinction of two competing species subject to contest competition and Allee effects. 


5.1. Two dimensional strong Allee effects of symmetric competition models 

Adopting the definition of Allee effects for single species, we are able to define Allee effects of a 
symmetric two species competition model that can be represented as 

x t+ i = F(x t ,y t ) = x t f(xt,yt) (18) 

yt +1 = F(y t ,x t ) = ytf(yt,x t ). (19) 

where < 0. Then we say two dimensional system (18)-(19) has two dimensional strong Allee 

effects if both f(x, 0) and f(x,x) satisfy Condition Al, A2, A3. Model (9)-(9) subject to inequalities 
(11) provides an example of two competing species suffering from contest competition and two dimen¬ 
sional strong Allee effects. The ecological examples when two competing species suffering from scramble 
competition and two dimensional strong Allee effects are the following two examples: 
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Basins of Attractions 



Initial conditions of species x 


Figure 14: Basins of attractions of (14)-(15) when r± = r 2 — 2.5, b\ = 62 = .05 and 5% = 5^ = 5: The black region is 
the basins of attractions of ( 0 , 2 / 2 ); The cyan region is the basins of attractions of (#2,0); The blue region is the basins of 
attractions of the interior attractor which is an interior equilibrium and the red region is the basins of attractions of (0, 0). 


Model I: Scramble competition with strong Allee effects due to mating limitations 


Xt+l 


Vt+i 


x t e r{1 ~ Xt) - aVt 


bx t 

1 + bx t 


y te r (i-m)- a x t 


byt 

1 + by t 


( 20 ) 

( 21 ) 


where r represent the intrinsic growth rates as well as intra-specific competition coefficient; b represent 
the product of an individual’s searching efficiency and the carrying capacities and a represents the 
inter-specific competition coefficient. If 


1 


- > 

r + a r + a 


1 

b’ 


then Model I has two dimensional strong Allee effects. 


Model II: Scramble competition with strong Allee effects due to predation saturations 


x t +1 

yt+i 


x t e r{1 - Xt) - ayt -^i 

y t e r{1 ~ Vt) ~ aXt ~ 


( 22 ) 

(23) 


where r represents the intrinsic growth rates as well as intra-specific competition coefficient; m represents 
predation intensities; b represent the product of the proportional to the handling time and the carrying 
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Nullclines 



Figure 15: Nucllines of (14)-(15) when r± = r 2 = 2.5, b\ = 62 = .01 and <$3 = £4 = 5. 


capacities and a represents the inter-specific competition coefficient. If 

(br + r + a) 2 


r < m < 


4b(a + r) 

then Model I has two dimensional strong Allee effects. 


and b > 1 + a/r, 


Kang and Yakubu (2011) and Kang (2013) have studied the detailed dynamics of Model I & II 
introduced above. By comparing our results to their study on a two competing species model where 
each species is subject to scramble competition and strong Allee effects (i.e., Model I &: II), we are 
able to obtain some generic dynamical features of a two-species competition model subject to strong 
Allee effects and the different dynamical outcomes due to different types of competition introduced in 
the models as follows: 


Generic dynamic features of scramble v.s. contest competition models with strong Allee 
effects: 

1. Both contest and scramble competition models with strong Allee effects can have either tri-stability 
(i.e., three boundary attractors Eq, E Xk ^, Eo yki ) or four attractors (i.e., three boundary attractors 
plus an interior attractor) depending on certain conditions. 

2. The coexistence of two competing species is possible only if the system has at least four interior 
equilibria. 

3. Both models can exhibit two dimensional strong Allee effects in the symmetric case. 
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Figure 16: Basins of attractions of (14)-(15) when r\ = r '2 = 2.5, b\ = h> = .01 and 83 = 84 = 5: The black region is 
the basins of attractions of (0, y 2 ); The cyan region is the basins of attractions of (x 2 ,0): The blue region is the basins of 
attractions of the interior attractor which is an interior equilibrium and the red region is the basins of attractions of (0, 0). 


4. The basins of attractions of Eq are unbounded, which suggests that two competing species are not 
able to coexist if their population densities are too low or too high. 

Different dynamical outcomes of scramble v.s. contest competition models with strong 
Allee effects: 

1. Our generalized Beverton-Holt two species (representing contest competition) model always has 
three boundary attractors E 0 , E Xki0 , E 0yki . While two-species Ricker models (representing scram¬ 
ble competition) subject to strong Allee effects can have either one (i.e., extinction of two species) 
or three boundary attractors where the extinction is caused by strong Allee effects (i.e., each species 
has essential extinction at its single state). 

2. Scramble competition models with strong Allee effects can promote the coexistence of two com¬ 
peting species at their high densities under certain condition. While this is not the case for contest 
competition models with strong Allee effects. 

3. Scramble competition models can have complicated structures of attractors (e.g., chaos, essential 
extinctions ) while the contest competition model may have relative simple dynamics (e.g., from 
performed simulations, we only observed locally asymptotically stable equilibriua as attractors). 

4. The basins of attractions of attractors that single species persists are bounded and may consist of 
different distinct parts for scramble competition models (see black and cyan areas in Figure 3&4 in 
Kang 2013) while the basins of attractions of attractors that single species persists are unbounded 
and connected for contest competition models (see black and cyan areas in Figure 12, 14 and 16). 
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The discussion above regarding contest v.s. scramble competition models with strong Allee effects 
are based on our study on the generalized Beverton-Holt competition model (7)-(8) and Ricker-type 
competition models with Allee effects studied in Kang and Yakubu (2011) and Kang (2013). The 
comparison may provide us useful insights of general two species competition models where each species 
has strong Allee effects in its single state. However, in order to obtain more general results, we need 
to classify the different types of intra- and inter-specific competition rigorously and study how their 
combinations with Allee effects lead to different dynamical outcomes. For example, the following example 
can be considered as a combination of contest inter-specific competition and scramble intra-specific 
competition subject to Allee effects: 


_ x t e ri(1 ^^(xt) 
t+1 1+ClVt 

_ yte r ^- y A p (y t ) 

y*+ X l+c 2 ;r t 

where I l ,i = x,y represents positive density dependence that can lead to Allee effects and Ci,i = 1,2 
are relatively inter-specific competition coefficients. This could be our future research project. 

5.2. Conclusion 

In this article, we propose a generalized Beverton-Holt competition model to investigate how the 
interactions of Allee effects and contest competition affect the species’s establishment and extinction. 
More precisely, we would like to answer the following questions addressed in the introduction: 1. What 
are the population dynamical outcomes of a two competing species when each of them suffers from both 
strong Allee effects and contest competition? 2. How does contest competition combined with strong 
Allee effects affect species’ persistence and extinction by comparing to models subject to strong Allee 
effects and scramble competition? 3. What are the generic dynamical features of population models 
subject to strong Allee effects and competition? Here, we summarize our main results from the study 
as follows: 

1. Theorem 3.1 and Theorem 3.2 reveal that a two species contest competition model with Allee 
effects has equilibrium dynamics, under the conditions that the intensity of intra-specific competi¬ 
tion intensity either exceeds that of inter-specific competition or is less than that of inter-specific 
competition but the maximum intrinsic growth rate of one species is not too large. This partially 
answers the first question. 

2. Theorem 3.3 and Theorem 3.4 provide us an approximations of basins of attractions of boundary 
attractors (i.e., initial conditions that lead to the extinction of one or both species). This indicates 
that the inter-specific competition can make species prone to extinction if it suffers from both strong 
Allee effects and contest competition, and the situation is worsen if species’ relative competition 
degree is less than 1. In addition, Theorem 3.5, 3.6 combined with numerical simulations suggest 
that decreasing the values of the inter-specific competition coefficient and degree may promote the 
coexistence of two species. These results may answer the second question. 

3. The third question can be partially answered by comparing with dynamics of two species scramble 
competition models subject to strong Allee effects in Kang (2013) to our Theorem 3.5 combined 
with numerical simulations. The comparison suggests that (i) Both contest and scramble compe¬ 
tition models can have either three boundary attractors or four attractor with only one interior 
attractor, however, scramble competition model may have only one attractor of the extinction of 
both species due to strong Allee effects; (ii) Scramble competition models represented by Ricker- 
type models may have more complicated structure of interior attractor than contest competition 
ones; and (iii) Scramble competition may be more likely to promote the coexistence of two species 
at low or high densities under certain conditions. 
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Our study combined with the results from Kang and Yakubo (2011) as well as Kang (2013) provide us 
useful insights on how Allee effects interacting with different forms of competition have an impact on 
the persistence and extinction of species. One direct application of these results is that we may be able 
to save rare species from essential extinction due to strong Allee effects by introducing another similar 
competing species under proper conditions. Thus, our study may help us understand the diversity of 
nature and make better policy for conservation program, especially for rare species that suffers from Allee 
effects. In the future, it will be interesting to investigate the following questions: 1. What are population 
dynamics of two competing species when each of them or both have different form of competition in 
their different life stages as well as Allee effects? 2. What are necessary conditions for a general two 
species population model to generate two dimensional strong Allee effects? 3. How can we predict more 
precise initial conditions that lead to the extinction and the coexistence of two species? 
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Appendix 

Proof of Proposition 2.1 

Proof. First, it is easy to see that Model (2) is positively invariant and bounded in M + . Notice that 

H( 0) = 0, 'AM = > 0 and Inn H(u) = r > 0, 

du (1 + U°y u->oo 

we can apply Proposition 2.1 of Kang and Dieter (2011) to obtain the following conclusion: Assume 
that u = H{u ) has n + 1 consecutive, distinct and non-degenerate solutions u l ,i = 0,1, ...,n with the 
following property 

0 = u° < u 1 < ■ ■ ■ < u n . 

then the even u l ,i > 2 is locally asymptotical stable with (u l_1 , u* +1 ) as its basin of attraction when 
u° = 0 is locally asymptotical stable (\.e. dH d ™’ | m < 1) and the odd u l ,i > 1 is locally asymptotical 
stable with (u 1-1 , tt I+1 ) as its basin of attraction when u° = 0 is unstable (i.e., dH d ^ | u > 1). 

According dH d ^ = , we can conclude that u° = 0 is unstable if <5 < 1 and u° = 0 is locally 

asymptotical stable if <5 > 1. Therefore, we have the following statement: i) If S £ (0,1) and Model (2) 
has two non-negative equilibria 0, K, then 0 is an unstable and and K is a stable positive equilibrium 
that attracts all points in (0,oo). ii) If <5 > 1 and the only equilibrium of Model (2) is 0, then 0 is 
globally stable, iii) If S > 1 and Model (2) has three non-negative distinct equilibria: 0, A and I\ such 
that 0 < A < K , then 0 equilibrium is locally asymptotical stable with the basin of attraction [0, A); A 
is unstable (repellor), while K is locally asymptotical stable with the basin of attraction (A, oo). 

Now we only need to derive the conditions when (2) has only one, or two or three equilibria. As we 
know, 0 is always a equilibrium of (2). If u* is a nontrivial positive equilibrium of (2), then we have 
h(u*) = 1, i.e., u* is a positive root of F(u) = u 5 — ra 5-1 + 1. If 6 < 1, then we have 

dF 

lim F(u) = —oo, lim F(u) = oo and —— = Su 5-1 — r(6 — 1 )u s ~ 2 > 0. 

u->0+ u-t oo du 
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This implies that F(u) = u has a unique positive solution u* and 


rS(u*) 


*\<5—1 


dH | 

d U 'u =“* (1 + (u*) s ) 2 1 + (u*) 5 


<5 < 1 . 


Therefore, if <5 £ (0,1), Model (2) has only two non-negative equilibria 0, A' where 0 is an unstable and 
and K is a stable positive equilibrium that attracts all points in (0,oo). 

If 5 > 1, then ^ = Su 5-1 — r(5 — 1 )u s ~ 2 = u s ~ 2 [Su — r(S — 1)] = 0 when u = u c = . Notice 

that 

dF, „ dF, 


F( 0) = 1 > 0, 


du ' u<u c 

thus F(u) =0 (i.e., h(u) = 1) have no positive root if 


< 0 and 


du u> “ c 


> 0 , 


5-1 r 


TVb 

F(u c ) > 0 => h(u c ) < 1 => c g 

1 +u° 


/ r(5-l) V 

—T 5 -< 1 => r < r crit = 6(6 - 1)?" 

+ (^A) 


and F(u) = 0 have two non-negative distinct equilibria 0 < A < u c < K if 

5-1 

> 1 => r > r crit = 6(6 - l)? -1 . 


,5-i r 


F(u c ) < 0 => h(u c ) > 1 => c , 

1 +u° 



Therefore the statements of Proposition 2.1 hold. 


□ 


Proof of Theorem 2.1 

Proof. First, we can easily verify that Model (6) is positively invariant and bounded in R + . If u* is a 
positive equilibrium of (6), then it is a positive root of F(u) = hu d + u s — rzi a_1 + 1. Since 

F'(u) = bdu d ~ l + 6u S ~ 1 - r(6 - l)^" 2 = u s - 2 [bdu d ~ 5+1 +6u- r(6 - 1)] = u s ~ 2 G(u) 
where G(u ) = bdu d ~ s+1 + 6u — r(6 — 1). If <5 < 1, then we have 

lim F(u) = — oo, F'(u) > 0, lim F'(u) = oo and lim F(u) = oo, 

u->0+ u-M3+ u-¥oo 

thus, F(u) has a unique positive root I\ > 0. 

Moreover, we have H( 0) =0 and H' (u) = ru ■ This indicates the following two cases: 

1. If 6 > d, then H' > 0, thus, Model (6) is a monotone system where its dynamics can be classified 
into the following two cases: 

(a) If <5 < 1, then Model (6) has two non-negative equilibria 0 and K where 0 is unstable and K 
is globally stable. 

(b) If 8 > 1 + d, then we have u = 0 is locally stable and there exists some unique positive u c such 
that G'(u) |u< u c < 0 and G'(u)\ u >u° > 0. Since lim u _ >0 + G(u) = oo and lim^oo G(u) = oo, 
therefore, G(u ) has exactly two positive roots if G(u c ) < 0. Also, we can conclude that 
G(u) > 0 if G(u c ) > 0. This implies that if G(u c ) < 0, then F(u) has exactly two positive 
distinct critical points 0 < u\ < u c < u 2 . Therefore, we can conclude that 1) F(u) has exact 
four distinct positive roots if F(u\) < 0 and F(u 2 ) < 0; 2) F(u) has exact two distinct positive 
roots if F(u\)F(u 2 ) < 0; and 3) F(u) has no positive root if F(u\) > 0 and F(u 2 ) > 0. In 
the case that F(u) has no positive root, we can conclude that u = 0 is globally stable. While 
if G(u c ) > 0, then F(u) < u for all u > 0, therefore, we have u = 0 is globally stable. 
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(c) Ifl + d><5>l, then we have u = 0 is locally stable and 

G( 0) = —r(S — 1), G'(u) > 0 and lim G(u) = oo, 

u—¥ oo 

thus, F'{u) has a unique positive point u c such that G(u c ) = 0 => F'{u c ) = 0 and 

F'{u)\ u<Uc > 0 and F'(u)\ u>Uc < 0. 

Therefore, if h[u c ) > 1, then Model (6) has as least two consecutive, distinct and non¬ 
degenerate equilibria A, K such that 0 < A < u c < K. 

2. If 6 < d, then there exists a unique critical point u c = ( r ( d ~ <5 ) )V rf suc h that 

H\u)\ u<Uc > 0, H'(u c ) = 0 and H'{u)\ u>Uo < 0. 

We consider the following two cases: 

(a) If 5 < 1, then Model (6) is a unique hump model that is similar to the Ricker’s map, thus, it 
has two non-negative equilibria 0 and K where 0 is unstable. 

(b) If 5 > 1, then H'( 0) = 0, thus u = 0 is locally asymptotically stable. If r > r cr it, then we 
have H(u c ) > u c . This implies that Model (6) has at least two distinct positive roots A. K 
such that 0 < A < u c < K. 

If both S and d are positive integers, then the number of interior equilibrium of (6) is the same as 
the number of positive roots of the following polynomial: 

F(x) = bx d + X s - rx 6 - 1 + 1. 

Descartes’ rule of signs (Smith and Latham 1954; Meserve 1982) states that if the terms of a single¬ 
variable polynomial with real coefficients are ordered by descending variable exponent, then the number 
of positive roots of the polynomial is either equal to the number of sign differences between consecutive 
nonzero coefficients, or is less than it by a multiple of 2. Therefore, the statement of Part 4 holds. 

□ 


Proof of Theorem 3.1 

Proof. Denote by T : the map defined by the right side of Model (7)-(8), then we have 


T 


(Ti(x,y),T 2 (x,y)) 


rix Sl 


1 + x Sl + biy Ss 


r 2 V S2 

1 + y &2 + b 2 x s 



(24) 


Thus, T is C 1 and Ti < r, t , i = 1, 2. The Jacobian matrix of T can be represented as follows: 


J 


rigi^i-Rl+bi/a) 

(l+j: 5 i +&iy 5 3 ) 2 

_ a,2b2S4y delta2 x 6 *- 1 

(l+y 5 2 +b 2 a; <5 4 


a 1 b 1 S 3 x s ^ y s 3 1 

(l+^i+bi y s 3 )' 2 

a2S2V 52 - 1 (l+b2X ° 4 ) 
(l+y 5 2 -|) 2 


(25) 


which gives 


det(J) 


rir 2 [<5i S 2 (x Sl 1 y 52 1 + b 2 x Sl+Si 1 y &2 1 +b\X Sl 1 y s ' 2+s 3 x ) + b 1 b 2 x 5l+5i 1 y S2+6 3 1 [5\8 2 

(1 + x Sl + biy Ss ) 2 (1 -f- x Sl + biy Ss ) 2 


- W] 


(26) 


According to Lemma 3.1, we know that the map T is positively invariant in IntM. 2 ^, i.e., T : Int Ry. —► 
InfR 2 + . Thus, for (x,y) £ Int R^_, we have det(J) > 0 (where equality occurs when Si 6 2 = 6364 ) and J 
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has a constant sign configuration as J = 


and it is -FT-strongly positive if S 1 S 2 > 8384 . 


. This indicates that (25) is K -positive in IntMf if 


Take u,v £ IntR 2 : such that u <k v, then by the mean value theorem and (25), we have 


Ti{v) — T\(u) = ^( Vl - Ul ) + ^(v 2 ~u 2 )>0 
T 2 {v) — T 2 (u) = ^{vi - ui) + %{V2 ~u 2 ) < 0 


where the partial derivatives are evaluated at a point on the line segment connecting u and v. This 
implies that T(u) «k T(v), i.e., the map T is strongly competitive. 

Choose any (xo,yo) £ T(IntR+), we have Xo < rq and yo < r 2 . Let (x,y) be a preimage of (xo,yo), 
then it satisfies the following equation: 


x 0 


yo 


nx Sl 

1 + x Sl + biy Ss 
T 2 V 62 

1 + y S2 + b 2 x s * 


(27) 

(28) 


This indicates that 


Xq 


yo 


l + x s 1+b tMl+hf 2 ly 3 /s 2 

1 v r 2 -yo > 

_ T2 y S2 _ 

l + y S 2 + b 2 ( X °^ 1 x yS3) ) S */ Sl 


fix) 

g(y)- 


(29) 

(30) 


Since <5 i<I 2 > 83 S 4 , we have 


fix) 

g'iv) 


ri(S 1 S2[x^+b 2 x^+ 5 i+b 1 x 5 H Voil r f 2 v l i) ) 63/S2 ]+bib 2 x 5 i + 5 *( m(1 +* 2 y l 4 Y 3/*2 [^ 1 g 2 -g 3 ^ 4 ]) ^ 


S 2 x(l+b 2 x s i)(l+x s i +b x ( 4 ) )^ 3/^2 )2 


r2(Si5 2 ly 62 +b iy 6 2 + 5 3+b 2 y 5 2 C° ( l+^ y 0 3 > ]+bi b 2 y s 2 +s 3 ( 31 )^ /l5 i [^i<5 2 — 6 3 ^ 4 ]) 


(31) 


^iy(l+biy s 3)(l+y s 2+b 2 (^±^Jl)S 4 /s l)2 


> 0. 


This implies that both / -1 (xo) and g~ 1 iyo) have only one point, i.e., ( xo,yo ) G T(IntM~i) has only one 
preimage. Thus, the map T is injective in IntMf^. 

It is easy to check that IntMf contains order intervals and is <x-convex (i.e., it contains the line 
segment jointing any two of its points that are A'-ordering). Therefore, we can apply Corollary 4.4 of 
Smith (1998) to obtain that {T n (xo,yo) = {x n ,y n )} is eventually componentwise monotone for every 
(xotVo) € IntMf, i.e., there exists a positive integer N such that either x n < x n +i for all n > N or 
x n +i < x n for all n > N and similarly for y n . Notice that 


T(IntR 2 + ) C T(R 2 + ) C [0,n] x [0,r 2 ], 


thus, {T n (xo, yo) = ix n ,y n )} converges to an equilibrium of Model (7)-(8) in [0,ri] x [0,r 2 ]. □ 


Proof of Theorem 3.2 

Proof. If 8 \ 8 2 < 8384 , then (26) vanishes at the following curves : 

LCf : xy = 0 

LC^r : (5 i<5 2 [1 + b 2 x Si + b\y S3 ] + bib 2 y S 3 x Si [ 8 2 8 i - $ 4 <S 3 ] = 0. 
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Thus, we have 


LC: y = 7 (a:) = 


( 5 i < 5 2 [1 + b 2 x 5 


1/^3 


or 


LCV X = 1 -\y) = 

which gives that 

Y 0*0 = - 


b 1 [b 2 x Si ( 848 3 - YYi) - <W 2 ] 

_<5i<5 2 [1 + hy Ss ] _ 

b 2 [hy S: > (8483 - 5 2 5 i) - <M 2 ] 

b 2 S iS 2 Slx Si ~ 1 


provided that x > 


(5i<5 2 


i/5 4 


, provided that y > 


5i 5 2 [1 + b 2 x 5 


^[6364 — ( SiJ 2 ] 


5 \S 2 

[^ 3^4 — <W 2 ] 


i/5 4 


i/5 3 


&i[6 2 ie' 54 ((>4<J3 — < 5 2 ( 5 i) — <M 2 ] 2 \61 \b 2 x Si (6463 — 8261) - 5 i 5 2 ] 

i/5 4 


i/s 3 —1 


< 0 


/ Ji \ 4 /°4 

provided that x > ( bo ~s 1 s 2 ] ) • I n addition, by straight forward calculations, we have 

j"(x) > 0 if (53 > 84 > 1. 


The curve 7 ( 2 ) separates IntMf+ into two connected unbounded components. The one unbounded 
component containing (0, e) x (0, e) for very small e is denoted by and it has det(J) > 0. While the 
other unbounded component is denoted by f2_ and it has det(J) < 0. 

Let T be the map describing Model (7)-(8). The image of 7 ( 2 ) under the map T, i.e., T(y), called 
the critical curves in (Mira et al 1996; Smith 1998; Kang 2011; Kang and Smith 2011), together with the 
portions of the coordinate axes which it cuts off, form the boundary of T (Int Mi). Because the curve 7 
is linearly order by «k (be-, 7 ^ 2 ) < 0 indicates that any two points in 7 are related in the sense of 
strictly K-ordering) and T preserves the strong ordering, T(y) is linearly ordered by «k and connects 
(0,r 2 ) to (ri,0). Let B denote the bounded component of IntM+ \ T(j). Then T(fi+) = T(fl_) = B, 
i.e., the map T folds Int over at 7 and maps both unbounded components onto B. 


If 


ri < 


SiS 2 


[<^3^4 ^ <5l<5 2 ] 


1/54 


or r 2 < 


5iS 2 


h[S 3 S4 — 5i<5 2 ] 


i/5 3 


then B is disjoint from 7 and is totally contained in f2+. This implies that det(T)B > 0. Note that the 
bounded component of Int-Mi \T( 7 ) contains the point without pre-image (e.g., choosing a point with 
x-coordinate larger than 77 and y-coordinate larger than r 2 ), then all the points in Int \ T{pf) has 
zero pre-image. This because that the cardinality of T~ 1 {w),w £ Int \ T(y) is finite and constant 
according to Lemma 3.4 of Chow and Hale (1982). Similarly, we can show that the cardinality of 
T~ 1 (w),w £ B is two where one is located in B and the other one is located in Int Mi. \ T( 7 ). This 
implies that 

T : B —> B is injective . 


Therefore, we apply Corollary 4.4 of Smith (1998) to obtain that {T n (x 0 , yo) = (x n , Vn)} is eventually 
componentwise monotone for every (xo,yo) £ Int R+. This indicates that every orbit of Model (7)-(8) 
with any initial condition in R 2 . converges to one of its equilibria. □ 


Proof of Theorem 3.3 

Proof. For any initial value (xo,yo) taken in O ex , we have Xq < Ai, thus we have 


xi = Ti(xo, yo) = -—jr < Ti(xo, 0). 

1 + Xq 1 + OlJ/ 0 


Since j£f\ x >o > 0, thus we have 


Xt +1 = T i(xo,yo) <T{{x 0 , 0 ). 
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According to Proposition 2.1, we have lim^oo Tf(xo, 0) = 0. Since Model (7)~(8) is positively invariant 
in by Lemma 3.1, thus, lirn^oo Xt = 0 holds for Model (7)-(8) if Condition HI is satisfied and its 
initial condition is taken in O ex . 

Similarly, we can use the same argument to show that lim^oo y t = 0 holds for Model (7)-(8) if 
Condition HI is satisfied and its initial condition is taken in O ey . Therefore, 

lim (x t , y t ) = (0, 0) = E 0 provided (x 0 , y 0 ) G O 0 . 


Based on our arguments above, we can conclude that for any initial value (xo,yo) taken in O x , we 
have limj-^oo y t = 0. In addition, we have 


x x = Ti(xo, yo) 


1 + Xg 1 + &i2/q 3 


> Ti(xq, A 2 ) 


nx S o 

1 + x* 1 +b 1 (A 2 y*' 


Then by induction, we have x t +x = Tl(xo,yo) > Tf(xo,A 2 ). Since Model (7)-(8) satisfies Condition H 2 
and xq > A“ x , then according to Corollary 2.1, we have 


lim Tl(x 0 ,A 2 ) = K J * 1 > AI 1 > Ax. 


This indicates that for Model (7)-(8), we have 


lim infer* > KV > A ? 1 > Ax- 

t—> OO 


Notice that lim^oo y t = 0, thus, for any e > 0, there exists a N large enough, such that for any t > N, 
we have 


x t +x > 


Si 

rxXt 1 

l + x^+e 


> A° r ' 


> Ax. 


Let e —> 0, then the limiting system of Model (7)-(8) becomes Ti(*,0). Since Condition HI is satisfied 
and the initial value is larger than Ax, thus according to Proposition 2.1, we have lim*-^ Xt = K\. 

Similarly, we can apply the same argument above to show the case when an initial value (xo,yo) is 
taken in O x . Therefore, the statement of Theorem 3.3 holds. □ 


Proof of Theorem 3.4 

Proof. Assume that Model (7)-(8) satisfies Condition HI. Let y = kx, then from (24) any point 
(x, kx) G IntM. 2 + maps to the following point 


T(x,kx) = (Tx(x,kx),T 2 (x,kx)) = 
If ^3 > <5i, then we have 


rxx Sl 


V 2 


( kx ) s2 


1 + x Sl + bi(kx) 53 ’ 1 + ( kx) s2 + b 2 x Si 
dTi(x,kx) r\x 5l—1 [<5i +bx(kx) S3 (Si —< 53 )] 


dx 


(1 + x Sl + bi(kx) S3 ) 2 


which indicates that 

dTi(x, kx) 
dx 


dTx(x,kx) 

x<x c ^ 0 cind |a:>x c ^ O 7 X c — 


Sx 


bxk S3 (5 3 -Sx) 


1/^3 


Let h{x) = a .i 1 +^(fc a .)i 3 = 1 + 6 lfc £ a; * 3 -«i > T i( x , kx )- Since K x ) > Tx{x,kx) and h'{x) < 0 provided 
(53 > <5i, thus we have 

f 1 
rx — Ax \ Ss ~ Sl 


T\{x,kx) < h(x) < Ai => x > 


bxk S3 
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This gives 


T\(x, kx) < Ai whenever x > 


ri - Ax 
bxk 


Define 


O l 


U«* , kx) G IntM. 2 ^ : x > 

k> 0 


f rx - Ax 
^ bxk s 3 



Then we can conclude that if Model (7)-(8) satisfies Condition HI and 6 3 > Sx, then for any initial value 
taken in O l ex , we have 

lim xt = 0 . 

£—>■ OO 

Similarly, we can conclude that if Model (7)-(8) satisfies Condition HI and ^4 > <5 2 , then for any initial 
value taken in Ol„ where 


o 


l 

ey 


1 J {{ky,y) G IntR 2 + : y > 

k> 0 


/ r 2 - A 2 
\ b 2 k Si 


1 



}, 


we have 


lim y t = 0 . 

£—>■00 


Let O q = O l ex n O l ey .Then based on the arguments above, we can conclude that if Model (7)-(8) satisfies 
Condition HI and <5 3 > > <5 2 , then for any initial value taken in O l 0 , we have 


lim ( x t ,y t ) = E 0 . 

t—> OO 

The argument above combined with the results of Theorem 3.3, we can conclude that the statement is 
hold. □ 


Proof of Theorem 3.5 

Proof. Since Model (7)-(8) satisfies Condition HI, thus if (x,y) is an interior equilibrium of (7)-(8), we 
have the the following two equations: 


1 = 

1 = 


Since 


thus we have, 


and 


where 


x+XT b ^ ^y=( rixH ~H XSl ~ 1 ) 3 =Fx{xlQ<Ax<x<Kx 
C'^ 2 -:/ 2 - 1 ) 1/54 = F 2 (y), 0 <A 2 <y< K 2 . 


riy A 

1+y d '2+6 2 a; d '4 


=> X = 


Fl(x) = 

F 2 (y) = 




1 / 83-1 


bid* 


y S 2 - 2 ( r2 y*2-l-yS 2 - 1 ) 

b 2 S 4 




[ri(5i - 1 ) - 8 xx) 

[r 2 (S 2 - 1) - S 2 y ], 


dFx , dFx dF 1 

7 \x<x c ^ O 5 7 tc=x c — 0 and x>a: c ^ O 5 

ax ax ax 


m dp 2 

I v<Vo > u 


dy 


dy 16 
rx{Sx - 1) 


= 0 and 


dF 2 


V>Vc 


dy 

r 2 (S 2 - 1) 


< 0 


and y c = 

Oi o 2 


Therefore, according to the geometry of y = F\ (x) and x = F 2 (y), we have follows: 
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• If Fi(x c ) < A 2 or F 2 (y c ) < A\, then Model (7)-(8) has no interior equilibrium. 

• If A 2 < F 1 (x c ) < K 2 , I< 1 < F 2 (y c ) or I< 2 < Fi(x c ), Ai < F 2 {y c ) < K 1 , then Model (7)-(8) has 
two interior equilibria. 

• If F\{x c ) > K 2 and F 2 (y c ) > K\. then Model (7)-(8) has four interior equilibria. 

□ 


Proof of Theorem 3 .6 

Proof. Let (x*,x*) be a symmetric interior equilibrium of Model (9)-(10), then we have 

r(sx *) 6 - 1 = 1 + {x*) s + b(x*) d . 

This indicates that the Jacobian matrix of Model (9)-(10) evaluated at the symmetric interior equilibrium 
(x*,x*) can be represented as 


J< 


( X*,X*) 


8(l+b(x*) d ) 

1 +(x*) 6 +b(x*) d 
bd,{x*) d 

l+(a:*) (5 +6(ai*) d 


bd(x*) d 

l+(x*) 6 -\-b(x*) d 

8(l+b(x*) d ) 

l+(x*) <5 +6(a:*) <i 


(32) 


The eigenvalues of (32) are 


AiOr*) 


6 + b(5 - d){x*) d 
1 + ( x*) s + b(x*) d 


and X 2 (x*) 


6 + b(S + d)(x*) d 
1 + ( x*) s + b(x*) d 


where Ai(a:*) denotes the stability of (x*,x*) on the invariant manifold f l y=x and the sign of A 2 (a:*) — 1 
denotes whether the eigenvector associated with A 2 is pointing towards or away from TL y=x . 

Therefore, the statement of Theorem 3.6 holds. □ 
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